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Multiwavelength observations suggest that clusters are reservoirs of vast amounts of relativistic electrons and 
positrons that are either injected into and accelerated directly in the intra-cluster medium, or produced as 
secondary pairs by cosmic ray ions scattering on ambient protons. In these possible scenarios gamma rays are 
produced either through electrons upscattering low-energy photons or by decay of neutral pions produced by 
hadronic interactions. In addition, the high mass-to-light ratios in clusters in combination with considerable 
Dark Matter (DM) overdensities makes them interesting targets for indirect DM searches with gamma rays. 
The resulting signals are different from known point sources or from diffuse emission and could possibly be 
detected with the Fermi-hAT. Both WIMP annihilation/decay spectra and cosmic ray induced emission are 
determined by universal parameters, which make a combined statistical likelihood analysis feasible. We present 
initial results of this analysis leading to limits on the DM annihilation cross section or decay time and on the 
hadron injection efficiency. 



1. Introduction 



Gamma rays from clusters of galaxies have not yet 
been detected^. On the other hand, from multiwave- 
length observations it is known that clusters have a 
very rich astrophysical phenomenology and observa- 
tions support the assumption of relativistic particle 
populations in the intra-cluster medium (ICM) . These 
particle populations can interact and produce high en- 
ergy gamma rays potentially being detectable by 
the Large Area Telescope (LAT) on board the Fermi 
satellite [Ij. 

Present day observations 0, |3] suggest the content 
of the universe to be made of about 25% of Dark Mat- 
ter (DM), whose nature however remains unknown. 
DM annihilation or decay into standard model par- 
ticles can lead to a gamma-ray flux, detectable by 
Fermi-L AT [|. The annihilation flux is proportional 
to the density squared, therefore galaxy clusters with 
substructure which can lead to large enhancements of 
the DM induced gamma-ray flux are promising targets 
for Dark Matter seaches. 

Both CR and DM-induced emission scenarios are 
different in spectral shape from any of the observed 
signals measured by the LAT, such as the emission 
from identified point sources or the galactic and ex- 
tragalactic diffuse emission. Assuming both the cos- 
mic ray (CR) and DM emission model being universal 
in at least one defining parameter in all clusters we 



^However, 
detected. 



galaxies who are members of clusters have been 



study, we can use a combined likelihood analysis to 
constrain the parameters that define those emission 
models. 

Both CR and DM induced gamma-ray emission 
need to be considered together in an attempt to quan- 
tify a possible DM or CR detection. In this paper, we 
present preliminary results where both components 
are considered separately. 

The outline of this proceedings is as follows. In Sec- 
tion 2 we discuss the cluster selection and their mod- 
eling with respect to the two emission models given 
by DM and CR, respectively. We present the analysis 
details and our combined likelihood method in Sec- 
tion 3. The results from the upper limit evaluations 
are described in Section 4. 



2. Cluster Selection and Modeling 

The clusters were selected based on the same crite- 
ria as in the first Fermi-LAT publication on DM in- 
duced gamma rays in galaxy clusters [Gj . We chose the 
brightest X-ray clusters from the extended HIFLUCS 
catalog 01 and removed clusters that were near the 
galactic plane (Norma and Ophiuchus) as well as clus- 
ters with detected active galaxies in gamma rays, such 
as Perseus with its central galaxy NGC 1275. In ad- 
dition we removed NGC 4636 since it exhibits a sig- 
nificant overlap in the region of interest with M49 and 
our method is at present not capable of treating that 
problem properly. 

It is worth mentioning that due to the lack of 
gamma-ray emission from any of these objects, we re- 
fer to the position in the sky which corresponds to 
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Figure 1: Predicted gamma-ray emission from neutral 
pion decay as in [ij]. We note that Coma is predicted to 
give the strongest signal. 



the NASA/IPAC Extragalactic Database locations of 
clusters as seen in other wavebands. 



2.2. Gamma-Ray Emission from DaxW 
Matter 

The gamma-ray flux from a pair-annhihilating DM 
particle with mass towimp can be written as: 
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• Bf X J A, 



(2) 



where {av) is the DM annhihilation cross section, 

and ■ Bf denotes the gamma-ray yield for the 

given final state with branching ratio B f . This parti- 
cle physics factor is multiplied with the astrophysical 
factor, J A that in case of annhihilating DM is the inte- 
gral over the line of sight and the DM density squared 
over a specified solid angle. For decaying DM we use 
the fact that the spectrum is roughly equivalent to 
the annhihilation spectrum of a particle with half the 
mass. Furthermore we note that the astrophysical fac- 
tor, Jd is only integrated over the DM density, not the 
squared one, in the case of decay. Eq.[2]then becomes: 



2.1. Gamma-Ray Emission from Cosmic 
Rays 
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where we denote the lifetime of the DM particle with 
r. 



Although there is no clear observational evidence 
yet for a relativistic proton population in clusters of 
galaxies, these objects are expected to contain sig- 
nificant populations of CR protons originating from 
different sources, such as structure formation shocks, 
radio galaxies, and supernovae driven galactic winds 
(see e.g. d, H, [IS [HI)- The cluster gamma-ray emis- 
sion is crucial in this respect as it potentially could 
provide unique and unambiguous evidence of a CR 
proton population in clusters through observing the 
pion bump in the gamma-ray spectrum. High res- 
olution simulations of galaxy clusters including CR 
physics show that these processes can be described 
by an approximately universal spectrum, where the 
gamma-ray flux from CR proton induced neutral pi- 
ous decaying can be parametrized as: 

^CB. ^ I a^rA{R)Ko^^ (E) (1) 

where AjrO^).^ (E) is the universal spectrum as given 
in [Tlj containing the maximum hadron injection effi- 
ciency, r], while A{R) is the individual normalization 
that varies with respect to the mass distribution as in- 
ferred from X-ray observations and its gas density 
profile [l2[- In fig. [T] we show the predicted gamma- 
ray emission from neutral pion decay for the clusters 
we chose to investigate. 



3. Analysis Details 

We carry out a binned analysis of 24 months of 
Fermi-LAT data. We select events corresponding to 
the DIFFUSE class which are most gamma-ray-like 
events and use the P6_V11 instrument response func- 
tion, provided in the software package Fermi Science 
Tools v9r22p0. The spectra are binned in 20 loga- 
rithmic bins from 200 MeV to 100 GeV. 

According to the model predictions for CR-induced 
gamma rays from galaxy clusters, the bulk of the emis- 
sion comes from the center of the cluster which justi- 
fies our choice of modeling them as pointlike sources. 
For the case of DM we assume the branching fraction 
to be 100 % to bb both for decay and annhihilation, 
which is shown to be a good approximation for a large 
class of possible DM models O [l3| • When modeling 
the clusters we model them as point sources using the 
same J-factors as in 6] for annhihilating DM and [lH| 
for decaying DM, without assuming any boost from 
substructure. 

We include all known point sources from the first 
year Fermi-LAT catalog [l^ that are within 15 de- 
grees from the location of each cluster using Fermi- 
LAT catalog fit results, except for point sources within 
5 degrees from the source, which are refitted. In ad- 
dition we model the diffuse and extragalactic gamma- 
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ray background according to the recommended Fermi- 
LAT templates^. 

When we inspect the spectra and the residuals of in- 
dividual sources with respect to the assumed model, 
we conclude that the measured 7-ray emission is mod- 
eled well by assuming only detected Fermi-LAT point 
sources and diffuse emission provided by the galactic 
and extragalactic diffuse models respectively. 



3.1. Likelihood Procedure 

In order to model each of the cluster regions given 
its own characteristics together with the underlying 
universal physics as assumed for emission from DM 
or CR respectively, we construct a combined likeli- 
hood function by multipylying every individual likeli- 
hood function and minimizing it with respect to the 
parameter of interest, e.g. the annhihilation cross sec- 
tion, the decay time or the hadron injection efficiency. 
The likelihood function can then be written as: 

L((o-u),mwiMp|o5s) 

= WLi{{(Tv),m'^iMi',Ci,hi\ohs) . 

i 

In Eq. |4] we use for simplicity the case of annhilating 
DM. For a given mass of a WIMP we can minimize this 
function and extract upper limits on the annhihilation 
cross section (crw). The 6, correspond to individual 
backgrounds that are treated as nuisance parameters 
in the likelihood evaluation, the Ci represent individual 
constants such as cluster masses. 

This method is implemented in the Fermi Science 
Tools as Composite2 and profiling over the likelihood 
function is achieved by means of MINOS which is part 
of the MINUIT package [13 • It should be stressed 
that we refer to the 95 % upper and lower limits as the 
ones that are defined as the points where the negative 
logarithm of likelihood function has changed its value 
by 2.71 with respect to its minimum. The method 
allows inclusion of the J factor uncertainties, see also 
for more details. In order to do so, the distri- 
butions of the J-factors need to be known. However, 
since the distribution of J-factors depends on under- 
lying assumptions such as whether or not DM sub- 
structure is included and the details of its assumed 
distribution and properties, we refrain from including 
J-factor uncertainties in our analysis. 



^For the galactic diffuse emission we use the template 
gll_iem_v02_P6_Vll_DIFFUSE.f it and for the isotropic gamma- 
ray background isotropic_iem_v02_P6_Vll_DIFFUSE.txt re- 
spectively. 
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Figure 2: 95 % upper limits on the Dark Matter 
annhihilation cross section for individual clusters (solid 
lines) and combined result (dotted) for bb final states. 
For comparison we also show the published values for the 
best cluster in the 11 months analysis (dashed) 01 . 




Figure 3: 95 % lower limits on Dark Matter decay time. 
Same as fig [2] 

4. Results 

We show the upper limits on the annhihilation cross 
section for our cluster sample and the combined re- 
sult in fig. [2] Below 100 GeV we exclude gamma rays 
from annhihilating DM with a annhihilation cross sec- 
tion above 5 x 10~^^cm'^s~^, while we note that our 
combined limit (dashed red line) is roughly a factor 
3 more constraining than the best individual cluster 
(Fornax). For the case of decaying DM we exclude 
life times below 10^*" s for WIMPs decaying into bb 
final states with masses above 100 GeV, as we show in 
fig [31 The adopted CR model assumes the maximum 
injection efficiency r; to be 50 %. We find that while 
only the case of Coma can constrain this parameter 
space, since the limit we find is below 50 %, that our 
data disfavors this value, see fig. ID Assuming a linear 
relation between the CR-induced gamma-ray flux and 
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Figure 4: 95 % upper limits on maximum injection 
efficiency following the hadronic cosmic ray model by 
[Ti]. The rightmost point indicates the combined 
constraint. 



rj, we can establish a combined upper limit on 77 of 
27 %, given the model characteristics by [111] , 



5. Summary and Conclusions 

Clusters of galaxies provide very interesting tar- 
gets for searches for gamma-ray emission from both 
hadronic CR interactions in the ICM and from Dark 
Matter annhihilation or decay. However, there is no 
confirmed observation of gamma rays from clusters so 
far. 

Individual fits of fermi-LAT data in the region 
of 5 clusters is compatible with identified Fermi- 
LAT point sources and galactic and extragalactic dif- 
fuse emission. From the non-detection of clusters in 
gamma rays, we can obtain upper limits on Dark Mat- 
ter annhihilation cross section or decay time for a 
generic WIMP with a bb final state. For the discussed 
emission model from hadronic cosmic ray interactions 
we are able to set upper limits on the hadron injection 
efficiency. 

To capture each region and its individual back- 
grounds we employ a combined likelihood analysis 
treating those individual backgrounds as nuisance pa- 
rameters. This approach is feasible as we assume all 
clusters to exhibit the same physical properties, which 
is undisputed for the DM case and can be assumed for 
the cosmic-ray emission to first approximation. Our 
combined limits on DM are a factor two better than 
previously published limits and our data disfavor CR 
hadrons being accelerated with an efficiency greater 
than 30 %. We note that one cannot lower the ac- 
celeration efficiency infinitely if one wants to explain 
radio halos with the hadronic CR model. 

It is however noteworthy that this analysis is only a 



proof-of-concept study to show the power of the com- 
bined likelihood analysis. We are currently updat- 
ing the analysis to increase the number of observed 
clusters, including more accurate spatial modeling, 
improved treatement of overlapping ROIs and explo- 
ration of other emission scenarios. We are also us- 
ing an updated Fermi-hAT event selection. A corre- 
sponding Fermi-LiAT publication is in preparation. 
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